Fast non-negative deconvolution for spike train inference from population calcium imaging

Calcium imaging for observing spiking activity from large populations of neurons are quickly gaining popularity. While the raw data are fluorescence movies, the underlying spike trains are of interest. This work presents a fast non-negative deconvolution filter to infer the approximately most likely spike train for each neuron, given the fluorescence observations. This algorithm outperforms optimal linear deconvolution (Wiener filtering) on both simulated and biological data. The performance gains come from restricting the inferred spike trains to be positive (using an interior-point method), unlike the Wiener filter. The algorithm is fast enough that even when imaging over 100 neurons, inference can be performed on the set of all observed traces faster than real-time. Performing optimal spatial filtering on the images further refines the estimates. Importantly, all the parameters required to perform the inference can be estimated using only the fluorescence data, obviating the need to perform joint electrophysiological and imaging calibration experiments.Comment: 22 pages, 10 figure

Agreement between two large pan-cancer CRISPR-Cas9 gene dependency data sets.

Genome-scale CRISPR-Cas9 viability screens performed in cancer cell lines provide a systematic approach to identify cancer dependencies and new therapeutic targets. As multiple large-scale screens become available, a formal assessment of the reproducibility of these experiments becomes necessary. We analyze data from recently published pan-cancer CRISPR-Cas9 screens performed at the Broad and Sanger Institutes. Despite significant differences in experimental protocols and reagents, we find that the screen results are highly concordant across multiple metrics with both common and specific dependencies jointly identified across the two studies. Furthermore, robust biomarkers of gene dependency found in one data set are recovered in the other. Through further analysis and replication experiments at each institute, we show that batch effects are driven principally by two key experimental parameters: the reagent library and the assay length. These results indicate that the Broad and Sanger CRISPR-Cas9 viability screens yield robust and reproducible findings

A high-resolution integrated map of copy number polymorphisms within and between breeds of the modern domesticated dog

<p>Abstract</p> <p>Background</p> <p>Structural variation contributes to the rich genetic and phenotypic diversity of the modern domestic dog, <it>Canis lupus familiaris</it>, although compared to other organisms, catalogs of canine copy number variants (CNVs) are poorly defined. To this end, we developed a customized high-density tiling array across the canine genome and used it to discover CNVs in nine genetically diverse dogs and a gray wolf.</p> <p>Results</p> <p>In total, we identified 403 CNVs that overlap 401 genes, which are enriched for defense/immunity, oxidoreductase, protease, receptor, signaling molecule and transporter genes. Furthermore, we performed detailed comparisons between CNVs located within versus outside of segmental duplications (SDs) and find that CNVs in SDs are enriched for gene content and complexity. Finally, we compiled all known dog CNV regions and genotyped them with a custom aCGH chip in 61 dogs from 12 diverse breeds. These data allowed us to perform the first population genetics analysis of canine structural variation and identify CNVs that potentially contribute to breed specific traits.</p> <p>Conclusions</p> <p>Our comprehensive analysis of canine CNVs will be an important resource in genetically dissecting canine phenotypic and behavioral variation.</p

Statistical Analysis of Molecular Signal Recording

A molecular device that records time-varying signals would enable new approaches in neuroscience. We have recently proposed such a device, termed a “molecular ticker tape”, in which an engineered DNA polymerase (DNAP) writes time-varying signals into DNA in the form of nucleotide misincorporation patterns. Here, we define a theoretical framework quantifying the expected capabilities of molecular ticker tapes as a function of experimental parameters. We present a decoding algorithm for estimating time-dependent input signals, and DNAP kinetic parameters, directly from misincorporation rates as determined by sequencing. We explore the requirements for accurate signal decoding, particularly the constraints on (1) the polymerase biochemical parameters, and (2) the amplitude, temporal resolution, and duration of the time-varying input signals. Our results suggest that molecular recording devices with kinetic properties similar to natural polymerases could be used to perform experiments in which neural activity is compared across several experimental conditions, and that devices engineered by combining favorable biochemical properties from multiple known polymerases could potentially measure faster phenomena such as slow synchronization of neuronal oscillations. Sophisticated engineering of DNAPs is likely required to achieve molecular recording of neuronal activity with single-spike temporal resolution over experimentally relevant timescales.United States. Defense Advanced Research Projects Agency. Living Foundries ProgramGoogle (Firm)New York Stem Cell Foundation. Robertson Neuroscience Investigator AwardNational Institutes of Health (U.S.) (EUREKA Award 1R01NS075421)National Institutes of Health (U.S.) (Transformative R01 1R01GM104948)National Institutes of Health (U.S.) (Single Cell Grant 1 R01 EY023173)National Institutes of Health (U.S.) (Grant 1R01DA029639)National Institutes of Health (U.S.) (Grant 1R01NS067199)National Science Foundation (U.S.) (CAREER Award CBET 1053233)National Science Foundation (U.S.) (Grant EFRI0835878)National Science Foundation (U.S.) (Grant DMS1042134)Paul G. Allen Family Foundation (Distinguished Investigator in Neuroscience Award

Inference of patient-specific pathway activities from multi-dimensional cancer genomics data using PARADIGM

Motivation: High-throughput data is providing a comprehensive view of the molecular changes in cancer tissues. New technologies allow for the simultaneous genome-wide assay of the state of genome copy number variation, gene expression, DNA methylation and epigenetics of tumor samples and cancer cell lines

Adaptations to Submarine Hydrothermal Environments Exemplified by the Genome of Nautilia profundicola

Submarine hydrothermal vents are model systems for the Archaean Earth environment, and some sites maintain conditions that may have favored the formation and evolution of cellular life. Vents are typified by rapid fluctuations in temperature and redox potential that impose a strong selective pressure on resident microbial communities. Nautilia profundicola strain Am-H is a moderately thermophilic, deeply-branching Epsilonproteobacterium found free-living at hydrothermal vents and is a member of the microbial mass on the dorsal surface of vent polychaete, Alvinella pompejana. Analysis of the 1.7-Mbp genome of N. profundicola uncovered adaptations to the vent environment—some unique and some shared with other Epsilonproteobacterial genomes. The major findings included: (1) a diverse suite of hydrogenases coupled to a relatively simple electron transport chain, (2) numerous stress response systems, (3) a novel predicted nitrate assimilation pathway with hydroxylamine as a key intermediate, and (4) a gene (rgy) encoding the hallmark protein for hyperthermophilic growth, reverse gyrase. Additional experiments indicated that expression of rgy in strain Am-H was induced over 100-fold with a 20°C increase above the optimal growth temperature of this bacterium and that closely related rgy genes are present and expressed in bacterial communities residing in geographically distinct thermophilic environments. N. profundicola, therefore, is a model Epsilonproteobacterium that contains all the genes necessary for life in the extreme conditions widely believed to reflect those in the Archaean biosphere—anaerobic, sulfur, H2- and CO2-rich, with fluctuating redox potentials and temperatures. In addition, reverse gyrase appears to be an important and common adaptation for mesophiles and moderate thermophiles that inhabit ecological niches characterized by rapid and frequent temperature fluctuations and, as such, can no longer be considered a unique feature of hyperthermophiles